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ABSTRACT 



Context. The diagnostics of stellar flaring coronal loops have been so far largely based on the analysis of the decay phase. 
Aims. We derive new diagnostics from the analysis of the rise and peak phase of stellar flares. 

Methods. We release the assumption of full equilibrium of the flaring loop at the flare peak, according to the frequently observed delay between 
the temperature and the density maximum. From scaling laws and hydrodynamic simulations we derive diagnostic formulas as a function of 
observable quantities and times. 

Results. We obtain a diagnostic toolset related to the rise phase, including the loop length, density and aspect ratio. We discuss the limitations 
of this approach and find that the assumption of loop equilibrium in the analysis of the decay leads to a moderate overestimate of the loop 
length. A few relevant applications to previously analyzed stellar flares are shown. 

Conclusions. The analysis of the flare rise and peak phase complements and completes the analysis of the decay phase. 
Key words. Stars: flare - X-rays: stars - Stars: coronae 



1. Introduction 



Solar and stellar coronal flares are impulsive events well de- 
tected in the soft X-ray band. They are typically explained as 
due to the sudden increase of temperature and emission mea- 
sure of the plasma confined in single or groups of magnetic 
tubes (loops), caused by strong heat pulses. 



The flare light curves in the soft X-rays are typically char- 
acterized by a fast rise phase followed by a slower decay (e.g. 
Haisch et al. 1983). In the decay phase, the plasma cooling 
is due to radiation emission and to thermal conduction to the 
chromosphere. Since the cooling times for a confined plasma 
depend on the characteristics of the confining structure, and in 
particular on its length, the analysis of the decay phase has been 
extensively used to diagnose the flaring loops, and in particular 
their size (see Reale 2002 for an extensive review). This is im- 
portant not only because stellar flaring regions are unresolved, 
but also because, more in general, this is one of the few tools 
to obtain detailed information on the geometry of the stellar 
coronae (e.g. Favata et al. 2000a), and of any other phenomena 
involving plasma confined in magnetic structures. 



To summarize, Serio et al. (1991) derived a thermodynamic 
time scale for the pure cooling of flaring plasma confined in 
single coronal loops: 



3.7 X 10" 



120- 



L9 



(1) 



where L (Lg) the loop half-length (in units of 10'' cm), To (Tqj) 
the loop maximum temperature (in units of 10^ K). In princi- 
ple, the loop length can be derived simply by inverting Eq.([Tll. 
Jakimiec et al. (1992) and Sylwester et al. (1993) analyzed ex- 
tensively the decay of solar coronal flares and pointed out that 
significant heating can be released even during the late phases 
of a flare. This can be diagnosed through the analysis of the 
path of the decay in a density-temperature diagram: sustained 
heating slows down the plasma cooling but much less the den- 
sity decrease. If this effect is not properly taken into account, 
the loop length can be significantly overestimated. Reale et al. 
(1997) use the scale time Eq.([T]i to derive a formula for loop 
length, corrected to include the effect of significant heating in 
the decay (see the Appendix for a review). This approach has 
been extensively applied to analyze observed stellar flares (see 
Reale 2003 for a review), with the exception of flares where 
the residual heating completely drives the flare decay over the 
plasma cooling. These flares are more appropriately described 
with models of two-ribbon flares (Kopp & Poletto 1984). 
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Eq.([Tl) and all the subsequent derivations lie on the assump- 
tion that the flare decay starts when the loop is at hydrostatic 
and energy equilibrium. Extensive modeling of solar and stel- 
lar flares has shown that the heat pulses that trigger the flare 
are impulsive, i.e. their duration is small with respect to the 
overall flare duration. One then wonders whether the flaring 
loop has had enough time to reach equilibrium before the heat 
pulse is switched off or not, and, if not, whether this might be 
important for the analysis of the decay. Jakimiec et al. (1992) 
showed that temperature and density begin to decrease simul- 
taneously if the heating lasts long enough to reach equilibrium. 
In many flares, it is instead observed that the temperature peaks 
(and therefore begins to decrease) measurably before the emis- 
sion measure, both in solar flares (e.g. Sylwester et al. 1993) 
and in stellar flares (van den Oord et al. 1988, van den Oord 
et al. 1989, Favata et al. 1999, Favata et al. 2000b, Stelzer et 
al. 2002). Cargill & KUmchuk (2004) pointed out that in tran- 
siently heated loops the cooling is initially dominated by ther- 
mal conduction and that the density begins to decay as soon as 
the radiative and conduction cooling times become equal. 

Most attention has been so far dedicated to the decay phase, 
which is the longest-lasting part of flares and therefore typi- 
cally off'ers more opportunities of time-resolved data analysis 
and of good photon statistics. Little and partial attention has 
been instead devoted to the initial phase of the flares. Hawley 
et al. (1995) study the rise phase jointly to the decay phase, in- 
cluding neither the heating in the decay, nor the delay between 
the temperature and the density peak. 

Flare observations from recent missions such as Chandra 
and XMM-Newton often detect flares in great detail and suc- 
ceed in resolving the rise phase (e.g. Giidel et al. 2004). On 
the other hand, in long flares it happens that the observations 
ends early in the decay phase, inhibiting the related analysis 
and making any kind of information derivable from the rise 
phase important (e.g. Giardino et al. 2004). 

In this work, we address specifically the rise and peak 
phase of flares, and investigate what diagnostics can be ex- 
tracted from its analysis. We will study the flare initial phases 
as a stand-alone analysis, and compare and cross-check with 
the analysis of the decay. We will also address possible ad- 
ditional diagnostics, e.g. the density and the loop aspect ratio 
cross-section radius over length - are flaring loops fat or slim 
or arcades of loops? - which cannot be constrained just from 
the decay analysis. In our derivation, we will maintain the as- 
sumption that the flare occurs in a single loop. This assumption 
is more realistic in the rise phase, when the impulsive heating 
typically involves a dominant loop structure while later resid- 
ual heating may be released in other similar adjacent loops (e.g. 
Aschwanden & Alexander 2001). We will instead release the 
assumption that the flaring loop evolves to a condition of equi- 
librium, and therefore also address the question of what is the 
efifect of releasing equilibrium conditions on diagnostic formu- 
lae, i.e. if the decay starts before the loop reaches equilibrium 
conditions, and what is the error from assuming equilibrium 
conditions. We will derive the loop length from the rise phase. 

In Section |2] we analise the flare evolution with a general 
outline, operative definitions, relevant results from modeling. 
In Section |3] diagnostic tools for the analysis of the flare rise 
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Fig. 1. Scheme of the evolution of temperature (T, thick solid 
line), X-ray emission (LC, thinner solid line) and density (n, 
dashed line) during a flare triggered in a coronal loop by a heat 
pulse. The flare evolution is divided into four phases (I, II, III, 
IV, see text for further details). 

and peak phase are derived. In Section|4] we discuss the results, 
the limitations of the analysis, the applications, with some spe- 
cific examples and in Section|5]we draw our conclusions. 

2. Flare analysis 

2.1. General flare evolution 

We consider a flare occurring in a single coronal loop. A flaring 
coronal loop can be modelled as a magnetic flux tube where 
the plasma is heated to flare temperatures by a transient heat 
pulse. The plasma confined in the loop can be described as a 
compressible fluid which moves and transports energy along 
the magnetic field lines (e.g. Priest 1984). 

We will suppose that: (i) the flare occurs in a semicircular 
loop with half-length L, initially in equilibrium conditions at 
much lower temperature and density than at flaring conditions; 
(ii) the flare is triggered by a heat pulse uniformly distributed 
in the loop; (iii) the heat pulse is a top-hat function in time; (iv) 
there is no heating in the decay; (v) the flaring loop is much 
smaller than the pressure scale height at the flare temperature. 

The plasma cooling is governed by the thermal conduction 
to the cool chromosphere and by radiation from optically thin 
conditions. The evolution of the confined plasma is well-known 
from observations and from modeling (e.g. Nagai 1980, Peres 
et al. 1982, Cheng et al. 1983, Nagai & EmsHe 1984, Fisher et 
al. 1985, MacNeice 1986, Betta et al. 2001) and in the follow- 
ing we summarize it into four phases, sketched in Fig.[T] which 
map on the path drawn in the density-temperature diagram of 
Fig.|2](see also Jakimiec et al. 1992). 

Phase I : from the start of the heat pulse to the temperature 
peak {heating). The heat pulse is triggered in the loop 
and the heat is efficiently conducted down to the much 
cooler and denser chromosphere. The temperature rapidly 
increases in the whole loop, with a time scale given by the 
conduction time in a low density plasma (see below). 
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Phase II : from the temperature peak to the end of the heat 
pulse (evaporation). The temperature settles to the max- 
imum value (Tq). The chromospheric plasma is strongly 
heated and expands upwards, filling the loop with much 
denser plasma. The evaporation is explosive at first, with 
a timescale given by the isothermal sound crossing time: 

L Lg 



Tsd 



25- 



(2) 



where ks is the Boltzmann constant, m is the average par- 
ticle mass. After the evaporation front has reached the loop 
apex, the loop continues to fill more gently. The time scale 
during this more gradual evaporation is dictated by the time 
taken by the cooling rate to balance the heat input rate (see 
Sec.O. 

Phase III : from the end of the heat pulse to the density peak 
{conductive cooling). When the heat pulse stops, the plasma 
immediately starts to cool due to the efficient thermal con- 
duction (e.g. Cargill & Klimchuk 2004), with a time scale: 



3nckBT()L^ IQ.SncksL?' 



112 



kT, 



5/2 



50 



p5/2 
0,7 



(3) 



where ("cjo) is the particle density (10'" cm"^) at the end 
of the heat pulse, /c = 9 x 10"^ (c.g.s. units) is the thermal 
conductivity. 

The heat stop time can then be generally traced as the 
time at which the temperature begins to decrease signif- 
icantly and monotonically. While the conduction cooling 
dominates, the plasma evaporation is still going on and the 
density increasing. The efficiency of radiation cooling in- 
creases as well. On the other hand, the efficiency of con- 
duction cooling decreases with the temperature. 
Phase IV : from the density peak afterwards (Radiative cool- 
ing). As soon as the radiation cooling time becomes equal 
to the conduction cooling time (Cargill & Khmchuk 2004), 
the density reaches its maximum, and the loop depletion 
starts, slow at first and then progressively increasing. The 
pressure begins to decrease inside the loop, and is no longer 
able to sustain the plasma. In this phase, radiation becomes 
the dominant cooling mechanism, with the following time 
scale: 



3kBTi^_ 



'ikBTt, 



9x10 



n3/2 

3 ^ M,l 



(4) 



nMP(T) bTl,nM nM.io 

where Tm (Tm,i) is the temperature at the time of the den- 
sity maximum (10^ K), hm (ha/.K)) the maximum density 
(10'° cm"-'), P(T) the plasma emissivity per unit emission 
measure, expressed as: 

P(T) = bT" 

with b - 1 .5 X 1 0" and a - -1/2, for consistency with the 
parameters of the loop scaling laws (Rosner et al. 1978). In 
this phase, the density and the temperature both decrease 
monotonically. The presence of significant residual heating 
could make the decay slower. This can be diagnosed from 
the analysis of the slope of the decay path in the density- 
temperature diagram (Sylwester et al. 1993, Reale et al. 
1997, see the Appendix). 





n 

Fig. 2. Scheme of the flare evolution of Fig. [T]in a density- 
temperature diagram (solid line). The four phases are labelled. 
The locus of the equilibrium loops is shown (dashed-dotted 
line), as well as the flare path with an extremely long heat pulse 
(dashed line). The corresponding decay path (marked with EQ) 
is the one typically considered by flare analysis based on the 
decay phase. 



2.2. Heat pulse duration 

The evolution outlined in Sec. |2.1| concerns a flare driven by a 
transient heat pulse. The analysis of stellar flares based on the 
decay phase typically assumes that the decay starts from equi- 
librium conditions, i.e. departing from the the locus of the equi- 
librium loops with a given length (hereafter QSS line, Jakimiec 
et al. 1992) in Fig.|2] The link between this assumption and the 
flare evolution outlined above is shown in Fig. |2] if the heat 
pulse lasts long enough, phase II extends to the right, and the 
flaring loop asymptotically reaches equilibrium conditions, i.e. 
the horizontal line approaches the QSS line. It is worth not- 
ing that, if the decay starts from equilibrium conditions. Phase 
III is no longer present, and Phase II links directly to Phase 
IV. Therefore, there is no delay between the beginning of the 
temperature decay and the beginning of the density decay: the 
temperature and the density start to decrease simultaneously. 
Also, the decay will be dominated by radiation cooling, except 
for its very beginning (Serio et al. 1991). 

The analysis of the rise and peak phase has to include the 
presence of Phase III, and the delay between the temperature 
peak and the density peak. This delay is often observed both in 
solar flares (e.g. Sylwester et al. 1993) and in stellar flares (e.g. 
van den Oord et al. 1988, 1989, Favata et al. 2000b, Maggio 
et al. 2000, Stelzer et al. 2002). The presence of this delay is a 
signature of a relatively short heat pulse, or, in other words, of 
a decay starting from non-equilibrium conditions. 



2.3. Hydrodynamic modeling 

We now use hydrodynamic simulations to analyze more in de- 
tail the evolution of the rise and peak phase of a flare (Phases 
I to III). The Palermo-Harvard code solves the time-dependent 
hydrodynamic equations for the plasma confined in a loop to 
describe the evolution of the density, temperature and veloc- 
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ity of the plasma along the loop (Peres et al. 1982, Betta et al. 
1997). 

As a representative example, we consider an initially quiet 
coronal loop with half-length L - 2 x 10"^ cm and tempera- 
ture of about 3 MK. A transient heat pulse is injected in it at 
time t=0. The flare heat pulse lasts for a time them, is uniformly 
distributed in the loop, and is as intense (9 erg cm"^ s"') as to 
bring the loop to a temperature of ~ 20 MK. 

Fig. [3] shows the evolution of the plasma temperature, den- 
sity and pressure at the loop apex for three different durations 
of the heat pulse (i.e. f/^o, = 90 s O.Stj, theai - 170 s ^ Ts and 
theat - 500 s a; St,) and for non-stopping heating. As long as 
the heat pulse is on, the simulation results overlap; they differ 
only for the decay phase, when the heating is off. The decay 
of the simulations with longer-lasting heat pulses begins later 
As already mentioned in Sec. 12.11 the temperature begins to 
decrease as soon as the heat pulse stops. The density and the 
pressure, instead, both continue to increase and they reach their 
maximum well later. We note that the shorter the heat pulse, 
the longer is the delay between the beginning of the tempera- 
ture decay and the density maximum. For the longest-lasting 
heat pulse, the density and pressure values are very close to 
the asymptotic equilibrium values, estimated from loop scaling 
laws (Rosner et al. 1978). We have checked that the evolution 
is self-similar for a longer loop L - 10"^ cm at the same tem- 
perature, with the evolution times scaling as Tj. 

In the density-temperature diagram (Fig.|4]i, the flare evolu- 
tion for the different heat pulse durations is well in agreement 
with that sketched in Fig.|2] For short-lasting pulses, phase IV 
starts as soon as the path crosses the QSS curve. For ^ 3ts, 
phase II ends and the decay (phase IV) starts both very close to 
the QSS curve, while phase III is practically absent (Jakimiec 
et al. 1992). 

3. Diagnostics of the rise and peak phase 

We now derive simple diagnostic tools for the analysis of the 
rise and decay phase of a flare, taking advantage of detailed 
numerical modeling. We first observe that the maximum pos- 
sible duration of the rise phase is the time taken by the loop 
to reach equilibrium conditions under the action of a constant 
(flare) heating. The simulation results in Fig. [3] show that the 
time taken by the plasma to reach equilibrium conditions is 
much longer than the sound crossing time (Eq.||2l), which rules 
the very initial plasma evaporation. This is also clear in Fig.|5] 
which shows the evolution of the pressure at the loop apex in a 
linear scale: after f = 50 s, the rate of pressure enhancement be- 
comes more gentle. As mentioned in Sec. 12. II in this phase the 
dynamics become much less important and the interplay be- 
tween cooling and heating processes becomes dominant. The 
relevant time scale is therefore that reported in Eq. ([U. 

After detailed analysis of extensive numerical modeling of 
decaying flaring loops, we have checked that the decay time 
from Eq. ([T]) should be computed with a correction factor: 

t', = (5) 

where (f) ^ 1.3 to better fit the decay time measured from nu- 
merical models. 




200 400 600 800 1000 
Time (s) 



10000 




Time (s) 

Fig. 3. Evolution of the plasma temperature {top panel), density 
(middle) and pressure (bottom) at the loop apex, as computed 
from a hydrodynamic model of a 20 MK flare in a loop with 
half-length L = 2 x 10^ cm, and with three different durations 
of the heat pulse (labelled in the top panel, in units of t,) and 
with non-stopping heating. The asymptotic density and pres- 
sure values are also shown (dotted horizontal lines). The time 
of the density maximum is marked (black spots). 
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Fig. 4. Evolution of the flare of Fig. [3] in the density- 
temperature diagram. The locus of the equilibrium loops de- 
rived from loop scaling laws (QSS curve, dashed line) is also 
shown. 
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Fig. 5. Pressure evolution of the flare of Fig. [3]in a linear scale 
to show that most of the rise phase can be reasonably described 
with a linear trend (dashed lines). 

Hydrodynamic simulations confirm that the time required 
to reach full equilibrium scales as the loop cooling time (Tj), 
and, as shown for instance in Fig. |5](see also Jakimiec et al. 
1992), the time to reach flare steady-state equilibrium is: 



3t, ^ 2.3t', 



(6) 



We have verified that Eq.© holds for loops of diff'erent lengths. 
Fig. |6] shows the ratio of the time required to reach 97% 
of the pressure equilibrium value to the time computed from 



0.5 - 



0.0 L 



1 10 
Loop holf-length (10* cm) 

Fig. 6. Time to reach pressure equilibrium obtained from hy- 
drodynamic simulations: the ratio of the time at which the pres- 
sure is 97% of the equilibrium value to the time obtained from 
Eq. (|6]l vs loop length. The size of the data points is propor- 
tional to the heating rate (maximum temperature of 10, 20 and 
30 MK). 

Eq. (|6]l for simulations of flaring loops with three different loop 
lengths, and three different heating rates appropriate to reach 
the temperature of 10, 20 and 30 MK, respectively. The agree- 
ment is within 10%. 

For t > teq, the density asymptoticafly approaches the equi- 
librium value: 



no - 



^ 



= 1.3 X 10" — 



la^kgL L 
where a - 1.4 x lO"' (c.g.s. units), or 



(7) 



"0,10 



13 



T- 

^ 0,7 



as obtained from the loop scaling laws (Rosner et al. 1978) and 
the plasma equation of state. 

If the heat pulse stops before the loop reaches equilibrium 
conditions, the flare maximum density is lower than the value at 
equilibrium. Fig. |5] shows that, after the initial impulsive evap- 
oration on a time scale given by Eq.©, the later progressive 
pressure growth can be approximated with a linear trend. Since 
the temperature is almost constant in this phase (Fig.|3]l, we can 
as well approximate that the density increases linearly for most 
of the time. We can then estimate the value of the maximum 
density at the loop apex as: 



"M ~ "0- 



(8) 



where tu is the time at which the density maximum occurs, 
which can be measured directly from observations. Since it is 
reasonable to assume that the volume of the flaring region does 
not change much, at least on the relatively short time scale of 
the rise phase, the time of the maximum emission measure is a 
good proxy for tu- 

Phase III ranges between the time at which the heat pulse 
ends and the time of the density maximum. Fig. [3] shows that 
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the start time of the temperature decay marks well the end of 
the heat pulse. In many stellar flares, the temperature evolution 
is not well resolved, and we may take the time of the tempera- 
ture maximum as indicative of the end of the heat pulse. 

The time of the density maximum is also the time at which 
the decay path crosses the locus of the equilibrium loops (QSS 
curve). As reported in Sec. |2?T] CargiU & Klimchuk (2004) 
remarked that, at this time, the radiative cooling time is ex- 
actly equal to the conductive cooling time. By equating the time 
scales in Eqs. ^ and Q, we can then derive the temperature 
Tm at which the flare maximum density occurs: 



= 9 X lO-\nML) 



1/2 



(9) 



or 



Tmj - 0.28(nM,ioi9) 



1/2 



If a value for L is already available to us and we are able to 
measure Tm, we can derive hm from Eq.®, in alternative to 
Eq.®. 

Since phase III is dominated by conductive cooling, we de- 
rive its duration, i.e. the time from the end of the heat pulse to 
the density maximum, as 

Afo-M ~ Tclni/r (10) 
where 

•A^ — 

Tm 

and Tf (Eq. 13)) is computed for an appropriate value of the den- 
sity «c- A good consistency with numerical simulations is ob- 
tained for = (2/3)nM- 

The collection of Eqs. (|6]l - (fTOl i provides a set of diagnostic 
tools for the analysis of the rise and peak phase of stellar flares. 
Eqs. (|6]l - ([8]l are related to the rise phase, the others to the peak 
phase, or, more precisely, phase III as defined in Sec. 12.11 We 
have checked that the equations provide values consistent with 
those obtained from accurate numerical modeling within few 
percent. 

By combining Eq.lfTOli and Eq.® we obtain: 



Afo- 



M 



1.21ni/r 



(11) 



which ranges between 0.2 and 0.8 for typical values of ifj (1.2 
-2). 

By combining Eqs.©, ([8]l and (|9]l, we derive a new expres- 
sion for the loop half length: 



L9 a! 3 lJ/^T^/^^tM,3 



(12) 



where fM,3 is Tm in units of lO"* s. For typical values of i// and 
To, we obtain Lg ~ 5 - 25fM,3 and Lg ~ lOfM.s may be taken 
for rough estimates. Therefore, we expect that flares occurring 
in loops with length of the order of 10'° cm show the peak of 
the emission measure about 1 ks after the beginning. 

By including Eqs. lfTTl i into Eq. (fT2l i. we derive another al- 
ternative expression for the loop half length: 



4,^ 



2.5 r^/2^ro-M,3 



where Afo-M,3 is AJq-m in units of 10-^ s. For typical values of 
tfr and To, we obtain Lg ~ 20 - lOOAfo-^.s and Lg ~ 50Afo-M.3 
may be taken for rough estimates. 

The loop length derived from application of Eqs.(fT2li and 
(fT3T l to the model simulations is correct within 10%. Eqs.(fT2]i 
and ( fT3] l allow us to estimate the length of the flaring loop from 
measuring characteristic time intervals of the flare rise phase 
and related temperatures. 

4. Discussion: implications and applications 

4.1. Consistency and limitations 

The aim of this work is to investigate the diagnostics of the 
rise and peak phase of coronal flares, to complement the well- 
established diagnostics of the decay phase (e.g. Reale et al. 
1997). Once derived the relevant diagnostic expressions, we 
first check whether there are effects on the analysis of the de- 
cay. In particular, we wonder whether assuming that the decay 
starts from loop equilibrium conditions - therefore ignoring the 
details of the "previous history" - leads to significant system- 
atic errors or not. In the decay analysis, one important param- 
eter is the temperature at the start of the decay. At equilibrium 
conditions, the decay starts at the temperature maximum. In 
Sec. 12.11 we have pointed out that, in an impulsive flare event, 
in which the loop does not reach equilibrium conditions, the 
density begins to decay later than the temperature. The proper 
decay phase begins at the density peak (i.e. the later time), and 
one should then use the temperature at the density peak, lower 
than the maximum temperature. The proper loop decay time 
then becomes: 



(14) 



If we use the maximum temperature Tq to estimate the loop 
length with expressions derived from Eq. (HJ (Reale et al. 
1997), instead of Tm, we then overestimate the loop length. 
However, since the dependence on the temperature is rather 
weak in Eq.® and the temperatures not very different, we ex- 
pect not a too large effect on previous results. The error can be 
easily estimated from the square root of the ratio of the maxi- 
mum temperature to the temperature at the density maximum. 
Since this ratio is typically of about 1.2-1.5 (see Sec l4.5l l, the 
amount of overestimate is no more than 15-20%. Furthermore, 
if we know this ratio we can correct for this effect. 

Fig.|7]shows how the decay time varies with the duration of 
the heat pulse. As expected, the decay time is invariably larger 
than the equilibrium decay time: the shorter the heat pulse, the 
longer is the decay time, with a maximum of about 1.7 t'^. 
However, in most of the range, the decay time is larger no more 
than 20% than the equilibrium decay time, and therefore the 
loop length overestimated by the same amount. 

In all the above analysis we have neglected the effect of the 
star gravity. This is reasonable both because the gravity have 
very little influence on the strong initial flare dynamics, and 
because, at the high temperatures of stellar flares, the pressure 
scale height: 



(13) h„~ 5000 
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Fig. 7. Time of the density maximum {tu, solid line), delay of 
the beginning of the density decay from that of the temperature 
decay (Afo-M , dashed line), and loop decay time (t^, dotted 
line) vs the duration of the heat pulse All times are in 
units of the equihbrium decay time t^. 

is larger than the height of the flaring loops (assuming they 
stand vertical on the stellar surface). For very long loops, com- 
parable to the stellar radius, the pressure scale height is even 
longer because the gravity decreases significantly at long dis- 
tances from the surface, so that the effect of gravity is expected 
to be even smaller. 

Finally, throughout our analysis we have assumed a sim- 
ple heating function: a top-hat in time, uniform distribution in 
space. Different heating functions may have some effects on the 
results (somewhat discussed in Sec. 14.5b . but except for very 
extreme cases, most of them should be smoothed out by the 
very efficient thermal conduction. 

4.2. Loop length from the rise phase 

Stellar flare observations do not always cover the whole flare 
duration; sometimes the rise phase is missing, sometimes most 
of the decay is not observed. The latter possibility can occur 
quite frequently, because the decay is the longest part of the 
flare and, for long flares, the observation may end too early in 
the decay. One result of the present work is to provide a set of 
useful formulae for events with missing observed decay. This 
could be obtained ultimately because most of the rise phase is 
ruled by the same processes which rule the decay phase, i.e. the 
energy losses. Eq. (fT2b yields the loop length, from the time of 
flie maximum density tM (the emission measure can be used as 
a proxy), the maximum temperature To and the temperature at 
the density maximum Tm- 

Alternatively, we may derive the loop length from Eq.(fT3ll 
even if we miss the flare start but we are able to measure the 
time interval between the temperature maximum and the den- 
sity (or emission measure) maximum. If both times are avail- 
able to us, we may choose the less uncertain one as the more 
reliable, and we may use both of them to check for consistency 
and to derive a more accurate length value from a weighted 



average. We should in fact consider that the determination of 
the relevant flare times can be aff'ected by significant uncertain- 
ties. The time of the density maximum is typically determined 
within a time bin where the spectral fitting is performed. This 
time bin can be quite large. Moreover, it is reckoned from the 
time of the flare start, which can in turn be not well-determined. 
Analogously, the uncertainty in the interval Afo-M comes from 
the width of both the time bins including the temperature maxi- 
mum and the density maximum. The strong dependence on the 
ratio if/ = To/Tm can also add to further uncertainties. 

Eqs. (fT2l i and (fTsT l provide the loop length even if only the 
rise and peak phase of the flare are observed. These expres- 
sions are alternative and independent of the expressions based 
on the flare decay and therefore of the presence of any signif- 
icant heating during the decay. They can provide therefore a 
further check on the loop length estimation, if both loop length 
values are available. There is a chance to obtain inconsistent 
results from the two approaches if the flare decay progressively 
involves other different coronal loops, or else if the heat pulse 
triggering the flare is not a top-hat function, instead, for in- 
stance, slowly rising. 

If no time-resolved spectral information is available, e.g. 
because of low photon statistics, the time of the light curve 
maximum may be used as a proxy of the time of the emission 
measure maximum. Since the former occurs a bit earlier than 
the latter, the loop length estimated from the former will be 
slightly underestimated. 

4.3. Loop Cross-section 

A typical output of the analysis of X-ray spectra with mod- 
erate resolution, such as those from CCD detectors (e.g. 
EPIC/XMM-Newton, ACIS/Chandra), is the emission measure 
associated with a fit thermal component. Time-resolved spec- 
tra during the flare can then provide us with the flare maximum 
emission measure. If we know also the maximum density, we 
can then derive the loop volume, and, knowing the loop length, 
also the loop cross-section area and radius. From Eq.® and/or 
Eq.(fT3b we can evaluate the maximum density at the loop top. 
Since the emission measure is integrated over most of the loop 
(the part emitting in the instrument band). For consistency, the 
volume should be computed using the density averaged of the 
emitting part of the loop. In realistic conditions of pressure 
equilibrium, this average density is higher than the density at 
the loop apex, because the temperature decreases downwards 
from the loop apex. A reasonable estimate of this average den- 
sity can be derived from the expressions linking the loop apex 
temperature to the temperature obtained from spectrum fitting 
(e.g. Favata et al. 2000b), under the reasonable assumption that 
the fit temperature is an average loop temperature: 

Tm^^ Tig (15) 

where the parameters ^ and 77 depend on the observation in- 
strument (e.g. ^ = 0.130 and 0.233 and 77 = 1.16 and 1.099 
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for EPIC/XMM-Newton, and MECS/BeppoS AX, respectively, 
see also the Appendix, Eq. |A.2| i. Then the average density is: 



Tm 



(16) 



Typically riayglriM ~ 1.5-2. The loop volume V, cross- 
section area A and radius R are then: 



EM 



avg 



2L 

A 
n 



(17) 
(18) 
(19) 



The final results of such formulae should be in any way 
taken with care, because highly indirect and therefore strongly 
affected by the error propagation. 

4.4. Single loop versus multi-loop 

It has been pointed out that large and long stellar flares likely 
involve entire coronal regions including multiple loop struc- 
tures. How can one reconcile this remark with the single loop 
approach followed in this work? Of course, since telescopes are 
not powerful enough to resolve the flaring structures, we can- 
not give conclusive answers. However, a few arguments can 
support our single loop approach. 

If multiple loop structures are involved in the flare, this fre- 
quently occurs in the late phases of a flare. The initial phases 
of an X-ray flare are instead quite localized and one can rea- 
sonably assume the presence of a single dominant loop, at least 
in the rise phase. This is observed in solar X-ray flares (e.g. 
Aschwanden & Alexander 2001) and supported by the accu- 
rate modeling of a well-observed stellar flare at least in the rise, 
peak and early decay (Reale et al. 2004). The clear evidence 
of a delay between the temperature and the density peak in 
many flares is consistent with a single loop model. Even in the 
later flare phases, a decay with no significant heating, i.e. with 
a steep path in the density-temperature diagram as sometimes 
observed even in very large flares (e.g. Favata et al. 2005) sug- 
gests strong plasma cooling confined in single loops. Arcades 
and two-ribbon flares are instead characterized by strong heat- 
ing which completely dominates the flare decay (e.g. Kopp & 
Poletto 1984) and/or by irregular light curves (Aschwanden 
& Alexander 2001, Reale et al. 200411 Recently, multi-thread 
models have been used to study solar flare evolution (Warren & 
Doschek 2005, Warren 2006). In this case, the hydrodynamics 
of the threads will be described by the results presented here. 



4.5. Sample applications 

As sample applications of the analysis described above, we 
have revisited three stellar flares already studied in the litera- 
ture: a flare on Algol observed with BeppoSAX on 30 August 



' Reale et al. (2004) showed that even arcades of equal loops can be 
described as single loops. 



1997 (Schmitt & Favata 1999, Favata & Schmitt 1999), a flare 
on AB Doradus observed with BeppoSAX on 29 November 
1997 (Maggio et al. 2000) and a flare on Proxima Centauri ob- 
served with XMM-Newton on 12 August 2001 (Gudel et al. 
2004, Reale et al. 2004). The first flare is quite long (~ 1 day) 
and shows an eclipse in the late decay phase. The first two 
flares are also hot (above 100 MK) and quite big, involving 
an emission measure above 10^"* cm"^. The last flare is cooler 
and on a smaller scale, but observed in great detail, thanks to 
the short distance of Proxima Centauri. It has been modelled 
accurately with detailed time-dependent hydrodynamic simu- 
lations (Reale et al. 2004), obtaining constraints even on the 
time and space distribution of the heating. These flares have 
been selected because for all of them the parameters related to 
the present analysis are all available with the associated uncer- 
tainties, and the uncertainties themselves are not excessively 
large. 

Table [T] shows the results obtained from Eqs.(|8Tl-(fT9]l for 
these three flares. The data in the first ten rows are derived 
from the original analysis of the references cited in the table. 
The first four are characteristic times, the light curve decay 
time (tlc), the time of the temperature maximum (fo), and of 
the density maximum (Im) reckoned since the beginning of the 
flare, and the delay between the temperature and the density 
maximum {dto-M)- There are then the flare maximum temper- 
ature (Tohs) and emission measure (EM) as derived from one- 
temperature fitting, and the ratio of the maximum temperature 
to the temperature at the density maximum. The loop maximum 
temperature (Tq), the half-length (L) and the thermodynamic 
cooling time (t^) are derived more indirectly according to em- 
pirical formulae (e.g. Reale et al. 1997, see the Appendix). 

The last ten rows show results obtained with the analysis 
presented in this work. The first four are densities: the loop 
equilibrium density («o) pertaining to the derived maximum 
temperature (Tq), the actual maximum density (««) derived in 
two different ways, and the maximum density averaged over 
the whole loop [tiavg), derived from the first maximum den- 
sity. Then there are the loop volume (V), cross-section area (A), 
and radius (r). Finally, we show other three values of the loop 
length: the first one is a refinement of the original length deriva- 
tion, using the temperature at the density maximum (Tm), in- 
stead of the maximum temperature (Tq). Since Tm < Tq, the 
new length is invariably smaller than the original, but not by 
large factors, as discussed in Sec l4.2l (within 20%), thus mostly 
confirming the results of the previous analyses. The other two 
length values are obtained from the analysis of the rise and peak 
phase, presented here (Eqs. {T71 and IfTSl ). 

The first two flares show a significant delay of the time of 
the density maximum from that of the temperature maximum 
- an indication that the heat pulse is relatively short as com- 
pared to the loop characteristic cooling time. Coherently, the 
ratio of the maximum temperature to the temperature at the 
density maximum is relatively large (e.g. Fig[3]). The flare on 
Proxima Centauri shows instead a relatively smaller delay and 
a coherently smaller temperature ratio. The loop length and the 
decay time are very large for the first flare (~ 10'^ cm and 21 
ks, respectively) and much smaller for the other two. 
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In all flares, the plasma does not reach equilibrium con- 
ditions, and the maximum density is well below the equilib- 
rium density. Although different, Eqs.® and (|9]l yield consis- 
tent density values - within the (quite large) uncertainties - for 
all flares, between a few 10'" cm"^ and a few lO" cm"^. The 
loop volume and cross-section parameters are affected by even 
larger uncertainties; they are derived from the combination of 
several quantities and suff'er from the error propagation. It is 
nevertheless confirmed the large loop aspect ratio for the AB 
Dor flare; the aspect ratio for the other flares is instead closer 
to typical solar coronal loops. 

The loop length values derived from Eqs. (fT2l l and (flJl i (re- 
ported in the last two rows) are generally consistent within 
the uncertainties with those derived from the analysis of the 
flare decay. They coherently yield quite smaller - although still 
marginally consistent - values for the AB Dor flare. We may 
take this as a vague indication that this flare involved progres- 
sively larger structures, coherently with the evidence of sig- 
nificant heating in the decay. We may also speculate that the 
opposite occurred in the Algol flare, i.e. an initial larger struc- 
ture, and later other smaller structures, consistent both with the 
significant heating in the decay and with the relatively smaller 
size obtained from the analysis of the eclipse. 

For the flare on Proxima Centauri, we obtain density val- 
ues smaller than those reported in Reale et al. (2004), which 
are, however, derived assuming quite a different heating depo- 
sition, i.e. concentrated at the loop footpoints. The loop aspect 
derived here is coherently larger than that reported for loop A 
in Reale et al. (2004), but if we consider the errors the results 
are almost compatible. All the loop lengths obtained for this 
flare are consistent with that constrained in Reale et al. (2004). 

We remark that in order to obtain consistent results, it is 
essential to take into careful account the uncertainties related 
to each step of the analysis. Table [T] shows that the analysis 
of the rise and peak phase of stellar flares provides valuable 
results and can therefore be usefully and extensively applied. 

5. Conclusions 

In this work, we derive useful diagnostic formulae for the anal- 
ysis of the rise and peak phase of stellar X-ray flares. The basic 
starting point is the realization that a flare is generally triggered 
by a short-lasting heat pulse, that the shorter the heat pulse the 
larger is the delay between the temperature maximum and the 
density maximum, and that the characteristic time scales in the 
rise phase also scale as the cooling times. We have then been 
able to derive useful expressions for the flaring loop density 
and length that can be applied if we can measure the flare rise 
times and a few significant temperatures. These expressions 
are generally independent of the analysis of the decay phase. 
Therefore, they can be used to complement and enrich the in- 
formation coming from the analysis of the decay phase, i.e. 
to check for consistency, to obtain constraints on the loop ge- 
ometry and even on the evolution of the flare morphology. Of 
course, they are particularly useful whenever the analysis of the 
decay phase cannot be performed, and we can equally derive an 
estimation of the loop length. Our analysis provides useful di- 
agnostics even when the data in the rise phase are limited. The 



information on the loop aspect represents a higher level of di- 
agnostics than that available only from the decay phase, and 
can therefore improve our knowledge on stellar flares and re- 
lated coronal structures. If a complete set of data is available, 
the complete analysis provides redundant information, and the 
opportunity of a cross-check; inconsistent results would not in- 
validate the present analysis, rather they would show that the 
related events are challenging because they cannot be well- 
described with our single-loop/simple-heat-pulse model. We 
look forward extensive application of this analysis to large sam- 
ples of stellar flares. The results presented here have application 
beyond flare loop evolution, such as the dynamical behavior in 
active region loops. 
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Table 1. Analysis of three stellar flaring loops. The parameters in the upper section are derived from the analysis in the references; 
those in the lower section from the equations presented in this work. 



Parameter 


Equation" 


Units 


AlgoI/BeppoSAX 


AB Dor/BeppoSAX 


Prox Cen/XMM 








(Favata & Schmitt 1999) 


(Maggio et al. 2000) 


(Reale et al. 2004) 


tlc 


dA.lb 


10' s 


50 ±5 


3.4 ±0.1 


4.3 ±0.1 


^0 




10 ^^ s 


13.5 ±4 


0.72 ± 18 


0.6 ±0.1 


tu 


® 


10 ^^ s 


34.5 ± 4 


1.12 ±0.09 


0.8 ±0.1 


A?0-M 




10' s 


21 ±6 


0.40 ±0.18 


0.19 ±0.09 


Tobs 




10' K 


14 ±2 


10 ±2 


2.6 ±0.1 


EM 




lO^* cm-3 


13 


5 


0.002 


Tq/Tm 


m 




1.4 ±0.2 


1.5 ±0.3 


1.18 ±0.07 


To 


llA.2t 


10' K 


20 ±3 


14 ± 3 


4.0 ± 0.2 




jAB 


10' cm 


800 ± 300 


49± 18 


13 ±4 


T? 


Q 


10' s 


21 ±8 


1.5 ±0.6 


0.8 ±0.2 


no 




10'" cm-' 


6.5 ±3 


50 ±20 


16 ±5 


Hm 


© 


10'" cm-' 


3.6 ±2 


12 ±6 


5±2 


flM 




10'" cm"' 


3.1 ± 1.2 


23± 12 


7±2 


flavg 


m 


10'" cm-' 


4.4 ± 1.8 


20 ±9 


9±2 


V 


Gl) 


10'" cm' 


7000 ± 5000 


330 ± 20 


0.25 ±0.1 


A 


GD 


10^" cm^ 


44 ±30 


34 ±24 


0.10 ±0.05 


r 


G2) 


10' cm 


37± 13 


33± 13 


1.8 ±0.5 


If 


iA.{\9\ 


10' cm 


670 ± 250 


40 ± 18 


12 ± 4 


L 


GD 


10' cm 


940 ± 300 


27 ± 13 


7±5 


L 


G2) 


10' cm 


1300 ± 700 


21 ± 14 


8±5 



" - where the parameter is used or evaluated; 

* - Time of the temperature maximum (or when the temperature begins to decrease); 

- The maximum temperature is used; 

- Using the original expression in Serio et al. (1991) with (p = I; 

- The temperature at the density maximum is used; 
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Appendix A: Review of the loop length from the 
flare decay 

From the basic work of Serio et al. (1991), and including the in- 
formation of the density-temperature diagram (Jakimiec et al. 1992, 
Sylwester et al. 1993) to diagnose residual heating, Reale et al. (1997) 
device a general expression of the loop length as a function of the ob- 
served decay time. The formula is essentially an inversion of Eq. G) 
with a factor (F(() > 1) which corrects for the presence of the heating: 



^min ^ ^ — Cm 



where <f is the slope of the decay in the log(n-T) (or equivalent ^EM- 
T) diagram and r^c is the e-folding time of the light curve decay (to 
1/10 from the maximum), and 



J- obs 



(A.2) 



Table A.l. Parameters for the determination of the flare loop 
length from the decay phase. 



Instrument 


Co 




9fl 








1 


ASCA/SIS 


61 


0.035 


0.59 


0.4 


1.7 


0.077 


1.19 


BeppoSAX/MECS 


0.68 


0.3 


0.7 


0.4 


1.7 


0.233 


1.099 


Chandra/AClS 


0.63 


0.32 


1.41 


0.32 


1.5 


0.068 


1.20 


EXOS AT/ME 


1.3 


0.4 


0.8 


0.4 


1.8 


0.195 


1.117 


ROSAT/PSPC 


3.67 


0.3 


1.61 


0.3 


1.8 


0.173 


1.163 


XMM/EPIC 


0.51 


0.35 


1.36 


0.35 


1.6 


0.130 


1.16 


GOES9 


1.02 


0.37 


0.36 


0.37 


1.7 


0.097 


1.163 


Yohkoh/SXT 


5.4 


0.25 


0.52 


0.3 


1.7 


(*) 


(*) 



(*) Yohkoh/SXT is supposed to resolve the flaring loop and to measure 
the temperature at the loop apex. 



and Tubs the maximum best-fit temperature derived from spectral fit- 
ting of the data. 

The correction factor is: 



(A.l) F{0 = 



(A.3) 



where c„, are parameters to be tuned for the instrument which 

observes the flare. 

Table IaTI shows the values of the parameters for some relevant 
solar and stellar instruments with moderate spectral capabilities. 



